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Abstract. - Fractal properties are usually characterized by means of various statistical tools 
which deal with spatial average quantities. Here we focus on the determination of fluctuations 
around the average counts and we develop a test for the study of galaxy distribution both in 
redshift and magnitude space. Fluctuations in the counts of galaxies, in a fractal distribution, 
are of the same order of the average number at all scales as a function of redshift and magnitude. 
We point out that the study of these kind of fluctuations can be a powerful test to understand 
the nature of galaxy clustering at very large scales. 



We consider the effect of real space long-range correlations in galaxy distribution on the 
statistical properties of counts both in redshift and in magnitude space. The existence of 
large scale structures (LSS) and voids in the distribution of galaxies up to several hundreds 
Megaparsecs is well known for twenty years [Q, 0]. The relationship among these structures 
on the statistics of galaxy distribution is usually inferred by applying the standard statistical 
analysis as introduced and developed by Peebles and coworkers || . Such an analysis assumes 
implicitly that the distribution is homogeneous at very small scale (Ao ~ 5 -j- 10h~ 1 Alpc). 
Therefore the system is characterized as having small fluctuations about a finite average 
density. If the galaxy distribution had a fractal nature the situation would be completely 
different. In this case the average density in finite samples is not a well defined quantity: it is 
strongly sample-dependent going to zero in the limit of an infinite volume. In such a situation it 
is not meaningful to study fluctuations around the average density extracted from sample data. 
The statistical properties of the distribution should then be studied in a completely different 
framework than the standard one. We have been working on this problem since some time [|] 
by following the original ideas of Pietronero || . The result is that galaxy structures are indeed 
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fractal up to tens of Megaparsecs g . Whether a crossover to homogeneity at a certain scale 
Ao, occurs or not (corresponding to the absence of voids of typical scale larger than Ao) is still a 
matter of debate M. At present, the problem is basically that the available redshift surveys do 
not sample scales larger than 50 -j- 100h~ 1 Mpc in a wide portion of the sky and in a complete 
way. In order to clarify this problem in the available galaxy catalogs, we propose here to study 
fluctuations around the average behavior of galaxy counts in redshift and magnitude space. As 
a matter of fact, in the counts of galaxies, large fluctuations around the average behavior have 
been reported by different authors both in redshift || and photometric surveys [ 
have been controversy as to whether they are due to real clustering or to incompleteness of the 
catalogs Q or to differences in the magnitude zero point of the various photometric surveys. 
It is, indeed, possible that discrepancies among these surveys are due mostly not to differences 
in photometric systems or in data reduction effects, but rather to real effects, i.e. large scale 
structures. We consider here this possibility and we present specific calculations about the 
expected value of these fluctuations in three cases: (i) a Poissonian homogeneous distribution, 
(ii) a purely fractal distribution, and (iii) a homogeneous distribution with scale-invariant (i.e. 
strongly correlated) fluctuations. 

Let n{r) = ^2iS(r — fl) be the microscopic density of a given set of galaxies. We assume 
hereafter that v(L, r) = <f)(L)n(r) is the density of galaxies in the volume element d 3 r around 
the point r, and with luminosity in the range [L,L + dL\. Note that by writing v(L,r) 
as a product of the space density for the luminosity function, the hypothesis that galaxy 
positions are independent of galaxy luminosity is implicitly assumed. Although it is known 
that there is a correlation between galaxy positions and luminosities jLq ], it has been tested 
that this is a reasonable assumption in the galaxy catalogs available so far |(|. The 
galaxy luminosity function usually assumed from observations is the so-called Schechter one: 
4>{L) — AL S exp(— L/ L„). The par ameters L* (luminosity cut-off) and 8 (power law index) can 
be determined experimentally |12| . The constant A is an overall normalizing factor such that 
the integral of <j>(L) over all luminosities is equal to one. In a volume limited (hereafter VL) 
sample W extracted from a given redshift survey, we may compute the number of galaxies as 
a function of distance of a fixed galaxy placed at the origin of coordinates. One can write the 
number-counts in a VL sample with faint luminosity limit at L — Lvl an d in a sphere C{R) 
of radius R in real (redshift) space as 

N(<R; >L VL )= f dL(f>(L) f d 3 rn(f). (1) 

JL VL JC(R) 

We may also consider the integrated galaxy counts as a function of apparent flux (or magni- 
tude), where / = L/4nr 2 is the apparent flux. In this case one simply have 

N(> f) = J™ dfn(f') = J™ dLc&L) J d 3 rn(r)Q - /) , (2) 

where Q(x) is the usual step function and the spatial integral is extended to the whole space. 
We study below average value and fluctuations of N(< r) (Eq.|l|) and N(> /) (Eq||) for an 
homogeneous, a fractal, and an homogeneous distribution with scale-invariant fluctuations. 

Before going on, some statistical definitions about averages and correlation functions have 
to be introduced kR. Consider a statistically homogeneous and isotropic particle density n{r) 
with or without correlations with a well defined average value uq. Statistical homogeneity and 
isotropy refer to the fact that any n-point statistical property of the system is a function only 
on the scalar relative distances between these n points. The existence of a well defined average 
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density means that 



C(R) 



(where ||C(i?)|| = 4ttR 3 /3 is the volume of the sphere C(R)) independently of the origin of 
coordinates. The scale A , such that | f c , R s d 3 r n(r)/\\C(R)\\ — n \ < n for R > A , is usually 
called homogeneity scale. If n(r) is extracted from a density ensemble, no is considered the 
same for each realization, i.e. it is a self-averaging quantity. Let (F) be the ensemble average 
of a quantity F related to n{r). If only one realization of n(r) is available, (F) can be evaluated 
as an average over all the different points (occupied or not) of the space taken as origin of 
the coordinates. The quantity (n(f{)n(f2) . . .n(fl)) dVidV2-.-dVi gives the average probability 
of finding I particles placed in the infinitesimal volumes dVi,dV 2 , ...,dVi respectively around 
^1)^2, ••■)'!■ For this reason (n(r\)n(r^ 2 )..m(fl)) is called complete /-point correlation function. 
Obviously (n(r)) = no, and in a single sample such that V 1 ^ Ao, it can be estimated by 
ny = N/V where N is the total number of particle in volume V. 

Let us analyze the auto-correlation properties of such a system. Due to the hypothesis of 
statistical homogeneity and isotropy, (n(ri)n(r2)) depends only on r\ 2 = \f[ — r$|. Moreover, 
(n(fl)n(f 2 )n(fi)) is only a function of r 12 = \f[ - f 2 \, r 23 = |ra - ^1 and r i3 = \f[ - r~i\. 
The reduced two-point and three correlation functions £(r) and C( r i2, ^23, ^13) are respectively 
defined by: 

(n(»!)n(T3)>=n§[l + £(r 12 )] (3) 
(n(fi)n(r 2 )n(ri)) = nl [1 + £(r 12 ) + ^(^23) + C( r i3) + C( r i2, ^23, ^13)] ■ 

In order to analyze observations from an occupied point it is necessary to define another kind 
of average: the conditional average (F) p . This is defined as an ensemble average with the 
condition that the origin of coordinates is an occupied point. When only one realization of 
n{f) is available, (F) can be evaluated averaging the quantity F over all the occupied points 
taken as origin of coordinates. (n(rl)n(r2)...n(r|)) p dVidVz—dVi is the average probability 
of finding I particles placed in the infinitesimal volumes dV\,dV2, ■•■,dVi respectively around 
r~i,r"2,...,n with the condition that the origin of coordinates is an occupied point. We call 
(n(fi)n(r~2)...n(fi)) conditional Z-point density. Applying the rules of conditional probability 
[ |l3| , one has: 

(n(f)) p = {n(6)n(r)/n (4) 
(n(fl)n(r" 2 )} p = (n(0)n(f{)n(fi)) /n . 

However, in general, the following convention is assumed in the definition of the conditional 
densities: the particle at the origin does not observe itself. Therefore (n(r)) p is defined only for 
r > 0, and (n(fl )n(f2)) p for r\, r-i > 0. In the following we use this convention as corresponding 
to the experimental data in galaxy catalogs. 

Imposing a partitioning on the space V (in the case V — > 00) with small cells of volume 
dV, let us consider the following Poissonian occupation process: 

„M = / rfV with probability n Q dV , . 

1 ' \ with probability 1 - n dV y ' 

where dV < 1/no and r is the center of dV. No correlation is supposed among different cells. 

1/3 

First of all one obtains (n(r)) = no and the homogeneity scale is given by Ao = 1/ ri ' . The lack 
of correlations implies (n(r\)n(f2)) = n\ if fx 7^ fh. and (n(fi)n(f2)) = no/dV if r\ =fh.. In 
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the limit dV — > one obtains £(?"i2) = S(f{ — f7)/nQ , Analogously, one can obtain the three 
point correlation functions: £( r i2j ^23, ^13) = 5{f{— f^) 6(1^—1%) /rig . The two previous relations 
say only that there is no correlation between different points. That is, the reduced correlation 
functions £ and £ have only the so called "diagonal" part. This diagonal part is present in 
the reduced correlation functions of any statistically homogeneous and isotropic distribution 
with correlations. For instance Jl4[ £(r) in general can be written as £(r) = 8(f) /uq + h(r) , 
where h(r) is the non-diagonal part which is meaningful only for r > 0. In general h(r) is a 
smooth function of r. As written before, in the definition of conditional densities, we exclude 
the contribution of the origin of coordinates. Consequently, we obtain: 

(11(f)) p = n Q (6) 
(n(fi)n(r 2 )) p = raj) [1 + S(f{ - fl)/n ] . 

By using Eq. (Q) and Eq. (||), one can compute both the average counts (N(< R; > Lvl))p and 
the relative fluctuation (AN 2 (< R; > L V l)) p = {N 2 (< R; > L V l)) p - (N(< R; > L V t))l , 
obtaining: 

(N(<R; >L VL ))p = r dL4>(L) (7) 

,, 2 , mN _ (N 2 (<R; >L VL )) p - (N(<R; >L VL )) 2 p 
((Wi?)>P= (N«R; >L VL ))l R 

This implies that the typical fluctuations between N(< R; > Lvl) seen by a single observer 

3 

and the average over all the possible observers goes to zero with R as R~? . In an analogous 
way, using Eq. (Q) and Eq. we can write: 

(N(> f)) p = AnaC(L*,5)f-% (8) 

(N 2 (>f)) P -(N(>f))l _ C 2 (L*,5) 1 f 
{ [I))p - (N(>f)) 2 p ~ C(U,5fn J ■ 

Then for / — ► the relative fluctuation of the /-counts from a single observer and the 
average behavior decreases as ft. The proportionality constants C(L*,$) and C2(L*,6) can 
be evaluated by performing explicitly the integrals implied by Eq. pi): 

C(L*,S) = l/(3V4^)i" +i r e (s + (9) 



C 2 (L* 1 5) = ATT(L sr ) 2S+ i J J™ dxdx'x s x' s e-( x+x '^ J™ dyy 2 e(^J^-y^ 6 




where T e (x) is Euler gamma-function. One can also consider the equivalent result in magnitude 
space. By definition § / = ilv{ ^' pc yi W 0A ( M *- m ) = KxlO-° Am . Then, from Eq.| one obtains 

(N(< m)) p ~ 10 a6m (10) 

(N 2 (<m)) p ~ (N(<m))l 

U) 2 m\\ =- - — - - - iiP^in- - 6 " 1 

{d {m))p ~ (N(< m)) 2 p W 

i.e. the normalized fluctuation exhibit an exponential decrease in m. 

Let us consider the fractal case. It is well known that fractal structures with the same 
fractal dimension may have very different morphological properties. How such properties 
can be defined and measured in a quantitative way is a very important question. Several 



Gabriclli & Sylos Labini, Fluctuations in galaxy counts 



5 



authors have introduced concepts like lacunarity, porosity, log-periodic corrections to scaling 
and others in order to go beyond the simple scaling analysis. Here we briefly review some 
useful concepts and definitions. For a fractal point distribution with dimension D < 3 the 
conditional one-point density (n(r)) p (which is hereafter called T(r)) has the following behavior 
[| (n(r)) = r(r) = Br D ~ 3 . Henceforth the average mass-length relation (hereafter MLR) 
from an occupied point is 

(N(< R)) p = (AnB)/D x R D , (11) 

The constant B is directly related to the lower cut-off of the distribution. Eq.(0) implies that 
the average density in a sphere of radius R around an occupied point scales as l/I$~ D . Hence it 
depends on the sample size R, the fractal is asymptotically empty and thus Ao — > oo. The MLR 
has a genuine power law behavior only when it is averaged over the ensemble, or over all the 
points of the structure if only one realization is available. The behavior of N(< R) from a single 
point presents scale invariant oscillations around the average given by Eq.([ll]). The presence 
of such oscillations is due to the fact that, at any scale R, the point distribution presents voids 
of radius of the same order as R. Usually this effect from a single point is described through 
a modulation term f(R) around the average behavior: N(< R) = (4irB)/D x R D x f(R) . 
\f(R)\ is a limited non-decreasing function of R. In the case of deterministic fractal f(R) can 
be written as [p| f(R) = exp (a sin(D/ log R + 4>) ) — 1 + a sin(D / log R + 4>) , with d«l. 
Therefore f(R) is an oscillating function with the wave length changing proportionally to the 
scale R. In a random fractal, f(R) has a more complex behavior, but keeps these qualitative 
features. By averaging over all the observers, f(R) is smoothed out: (f(R)) p = 1 . Note that 
the effect of f(R) is that fluctuations in N(< R) with respect to the average (N(< R)) p , are 
proportional to (N(< R)) p itself, rather than to its square root as in a poissonian distribution. 
Let us now consider the two-point conditional density (n(r~i)n(f2)} p . Blumenthal & Ball 
| jl7| have shown that statistical translation and rotational invariance, together with scale- 
invariance, lead to 

{n{r{)n{r 2 )) p ~ r(ri)r(r 2 )£(ri/rs, 9) , (12) 

where C(rx/r2, 0) — l+g(r\/r2, 9), with lixa ri / r2 ^o . r2 /r 1 -*o g( r i/ r 2,9) = 0, is called lacunarity 
function. 

By using Eq. (|l]), one obtains the average count in a VL sample: 

(N(<R; >L VL )) p = —R D x / d>{L)dL (13) 
u Jl vl 

{N 2 (<R; >L VL )) p - (N(<R; >L VL ))l ~ {N{< R; > L VL )) 2 p 

on large enough scales, when shot noise becomes negligible. Then, we expect that fluctuation 
from a single observer in spatial counts with respect to the average are of the same order of 
the average itself at any scale. Analogously, applying Eq. (||) to this case, one obtains 

(N(>f)) p =ABQ(L.,6)r% (14) 
(N 2 (> f)) P - (N(> f% = A 2 B 2 Q 2 (L^5)f- D ~ (N(> f)) 2 p , 

where 

C(L., S) = (l/(D(4^))Li +£ ^T e (s + (15) 

Q 2 {U, 5) = (L^ s + 2 + D JJ °° d Xl dx 2 xfx s 2 e-^+^ / °° d Vl y^ 1 
x f °° dy 2 y?- 1 J 4jr f 47r dn ld Q 2 g o) 9 - Vl ) Q y 2 ) 
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Eqs.(13)-([14|) show the persistent character of counts fluctuations in the fractal case. In terms 
of apparent magnitudes the normalized fluctuation becomes (8 2 (m)) p — Q,2{L*, 5)/[Q{L*, $)\ 2 . 
This implies that counts fluctuations around the average are constant with apparent magnitude 
(or with distance), in contrast with the previously considered poissonian case where they are 
exponentially (or power law) damped. A particular attention must be payed to the case D = 3 
in fact, from Eqs.(|il^-(|l4|), while the average counts formulas reduce to the Poissonian ones, 
fluctuations stay scale invariant in contrast with the Poissonian case. This is due to the 
assumed shape of the lacunarity function, which contains in itself the scale-invariant behavior 
of fluctuations. In order to obtain the real Poissonian-homogeneous case, it is not enough 
to consider the limit D — > 3, but it is necessary to substitute the relation between 2-point 
conditional density and I-point conditional density with {n(fi)n{r2)) p — {n{r\)) p (n(r2)) p at 
large enough \ f\ — r^l, and not large enough T\jr 2 which is a weaker condition. 

Let us now consider a mixed case in which the system is homogeneous (i.e. Ao is finite), 
but presents scale-invariant density fluctuations. This last event is in general described by 
the divergence of the correlation length r c , which is usually defined as the scale beyond 
which £(r) is exponentially damped. It measures up to which distance density fluctuations 
density are correlated. Note that while Ao refers to an one-point property of the system (the 
average density), r c refers to a two-points property (the density-density correlation) jl^, p"8|] . 
Therefore let us consider a system in which (n(r)) — no > and £(r) = [<5(r)]/no + f(r), with 
|/(r)| ~ r 1 for r ^S> \o , and < 7 < 3. The case 7 > 3 is not considered here. Indeed, since 
I J d 3 r£(r)| < +00, it is simple to show that fluctuation are at maximum Poisson-like, and in 
the case in which J d 3 r£(r) = 0, they can be also smaller. Moreover, this case corresponds to 
a scale-invariant fluctuation density field with "negative" dimension. 

We recall that in the definition of /-point conditional densities the contribution of the origin 
of coordinates is not considered. Therefore (n(r)) p = no[l + /(*")]• ^ n sucn a situation density 
fluctuations around the average have a scale- invariant nature. In fact a density field p(r), 
defined by p(r) — [n(f) — no], represents a kind of fractal field whose I-point conditional 
density is no£{r) ~ r~ 7 , and D = 3 — 7 represents the fractal dimension. For the same 
reason the quantity ((n(f{) — n ) (n(ra) — n )) = n 2 [£(ri 2 ) + C(ri, ?"2, ^12)] is analogous to 
the two-point conditional density of the fractal case. Therefore one can impose the equivalent 
of Eq.(|i~2|) to the present case: /(ra.2) + C( r i j r 27 ^12) = f{ r i)f( r 2)£( r i/ r 2:9) , where 9 is the 
angle between f[ and rj, and £(r 1 /r 2 ,6 | ) is defined by Eq.([l2[). Therefore, from Eqs. 
we can write 

(n(rl)n(rl)) p - (n(r"i)) p (n{r 2 )) p = n 6(fi - r 2 ) + nlf{ri)f(r 2 )g(ri/r 2 ,6) , (16) 

where the term in 5(fi — fi) is due to the diagonal part of £(ri2). As shown below, this 
term is important only for 7 > 3/2 (i.e. D < 3/2 if 7 > 0). In the case of a fractal point 
distribution this contribution was omitted because in that case it is always irrelevant. At this 
point we can evaluate counts fluctuations around the average both in a spatial (redshift) VL 
sample (for large R) and a flux limited one (for small /). First of all, we obtain (considering 
only the main contribution): (N(<R; >Lvl)) p ~ R 3 and (N(> /)) ~ . Moreover, 

if 7 < 3/2, we have {^ L {R)) P - R^ and (S 2 (f)) p - p , while if 7 > 3/2 the 
contribution from S(f{ — ri) in Eq.([l6|) dominates, and the same result as in the Poissonian 
case is obtained. Finally in terms of magnitude one find: (8 2 {m)) p ~ iQ-OAym f or ^ < 3^2 
and (5 2 (m)) p ~ lO -0 - 6 '™ for 7 > 3/2. Therefore also in the case of a homogeneous, 
but correlated point distribution, normalized counts fluctuations around the average must be 
always damped, no matter the large scale behavior of the correlation function. We have only 
a distinction of damping rate between the two cases 7 > 3/2 (which includes also the case 
of an exponentially tailed £(r)) and 7 < 3/2. In the first case, because of the fast decay of 
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£(r) at large scales, the behavior is the same as in the poissonian case without correlations; 
in the second case the damping of normalized fluctuations is slower as correlations become 
important. We conclude that the only case in which persistent and scale invariant normalized 
counts fluctuations can be observed, is the case of a fractal point distribution. Note that this 
behavior can be observed also in the case when the point distribution is homogeneous, but we 
are analyzing sample smaller than the homogeneity scale Ao- 

We have considered the average and fluctuations of the counts in the low redshift ap- 
proximation. The effect of the cosmological corrections (G: geometrical, K: k-corrections, E: 
luminosity evolution correction) are, in the range 0.1 < z < 1 linear with redshift. Such 
corrections, which are mostly model dependent, can change the slope of the average counts 
determined from a single point (i.e. The Earth). For example we have discussed the effect 
of K+G corrections in the determination of the slope of the counts as a function of distance 
from one point in [ p0[ . A similar effect can also occur in number counts as a function of 
magnitude in the equivalent regime (e.g. 15 m < B < 20 m ). However it is worth to note that 
the presence of possible persistent and scale-invariant fluctuations in the counts cannot be due 
to any smooth correction to data as cosmological and evolution effects, but they can be the 
outcome exclusively of the strongly correlated effects present in a fractal point distribution. 
In fact, smooth corrections to the average behavior cannot produce persistent scale- invariant 
fluctuations in the counts of the same order of the average itself. This property suggests an 
important experimental check of the possible fractal nature of the galaxy distribution from the 
behavior of number counts as a function of apparent magnitude with the advantage of using 
also two dimensional surveys, rather than only redshift ones. 
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